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Abstract. We investigate laser-induced nonsequential double ionization from 
aligned diatomic molecules, using the strong-field approximation in its length 
and velocity gauge formulations. Throughout, we consider that the first electron 
(^ ' dislodges the second by electron-impact ionization. Employing modified saddle- 

Q.' point equations, we single out the contributions of different scattering scenarios 

to the maxima and minima observed in the differential electron momentum 
distributions. We show that the quantum interference between the electron orbits 
starting and ending at a specific center Cj, and those starting at Cj and ending 
at a different center d, leads to the same maxima as minima as if all possible 
jrt ' scenarios are taken. There exist, however, quantitative differences as far as the 

• gauge choice is concerned. Indeed, while the velocity-gauge distributions obtained 
Jr^ , employing only the above-mentioned processes are practically identical to the 

overall distributions, their length-gauge counterparts exhibit an asymmetry in the 
(^ ' positive and negative momentum regions. This asymmetry is due to additional 

^~>' potential-energy shifts which are only present in the length-gauge formulation, 

<l ' and which, depending on the center, sink or increase the potential barrier through 

^ I which the first electron tunnels. In contrast, the interference between topologically 

similar scenarios leads at most to patterns whose positions, in momentum space, 
do not agree with the overall interference condition, neither in the length nor in 
C^ ' the velocity gauge. 
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Molecules in strong laser fields (/ ^ lO^"* — lO^^W/cm ) have been the object 
. , of intensive scrutiny in the past few years. In particular, high-order harmonic 

j^ ■ generation (HHG) and above-threshold ionization (ATI) have been employed to 

extract information about the structure of such systems with subfemtosecond precision 
[T] . This is made possible by the physical mechanisms behind such phenomena, namely 
the rescattering or recombination of the active electron with its parent molecule [5] ■ 
These processes take place within a fraction of a laser cycle. For a typical Titanium- 
Sapphire laser field, whose cycle is roughly T ~ 2.7 fs, this corresponds to hundreds 
of attoseconds [3]. In particular, the configuration of atoms with which the electron 
rescatters or recombines leads to patterns which are characteristic of the molecule. 
These patterns may be either due to the quantum interference of photoelectron 
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or harmonic emission in spatially separated centers j4j, or of the rescattering or 
recombination scenarios involving more than one center [51 [71 [HI [HI [13 [H] ■ 

Specifically for ATI and HHG in diatomic molecules, there exist analytic 
expressions which give the approximate energy positions of the interference minima 
and maxima due to the above-mentioned spatial separation [4]. This allows a 
physical interpretation in terms of a microscopic double-slit experiment. In particular 
in the framework of semi-analytic, S-Matrix approaches, such as the strong-field 
approximation (SFA), the HHG or ATI transition amplitudes are written as multiple 
integrals, with slowly- varying prefactors and a semiclassical action [1^. In this 
case, the double-slit interference appears as a prefactor, which depends on the 
symmetry of the highest occupied molecular orbital and on the internuclear distance 
|3[6l|7l[8l[9l[IIl[l3l[l4l[l5lll6l[18lll7]. 

There are, however, comparatively fewer studies of the influence of two-center 
scenarios, in which an electron is released in a center Cj in the molecule and rescatters 
or recombines with a different center C^, j ^ v [SJ [SJ [71 [HI [HI [IH]- In these references, 
the two-center processes have also been treated and discussed in quite different ways. 
For instance, in [71 [9l [10], the high-order harmonic double-slit prefactor has been 
exponentialized and incorporated in the action, while in high-order ATI the two-center 
processes led to transition amplitudes, which could not be grouped as to provide a 
common prefactor [S] . 

In previous work, we have shown that well-defined interference fringes may be also 
present in laser-induced nonsequential double ionization (NSDI) of diatomic molecules 
[H], for a small range of alignment angles. This is a consequence of the fact that 
NSDI is also described as a laser-induced rescattering process. More specifically, the 
first electron, freed at a time t' , is accelerated by the external laser field and, at 
a subsequent time t, collides inelastically with its parent molecule. In this collision, 
it gives part of the kinetic energy it acquired from the field to a second electron, 
which is then released. Quantum mechanically, transition amplitudes corresponding 
to NSDI at different centers in the molecule are expected to interfere. Even though, 
up to the present date, there is no experimental evidence of such fringes, recently, 
it has been observed that the shapes and the peaks of the NSDI differential electron 
momentum distributions in molecules depend on the symmetry of the highest occupied 
molecular orbital [20j , and on the molecular alignment angle with respect to the laser- 
field polarization [21]. 

In [19], we worked mainly within the SFA, and assumed that the second 
electron was dislodged by the simplest possible physical mechanism: electron-impact 
ionization. We have shown that, in this case, the two-center interference leads 
to minima and maxima parallel to the anti-diagonal P(i)|| -I- P(2)|| in the plane of 
the momentum components (p(i)||,P(2)||) parallel to the laser-field polarization. For 
parallel-aligned molecules, these fringes are sharpest. As the alignment angle increases, 
the contributions from the perpendicular momentum components start to blur these 
fringes, until, for perpendicular alignment, all structure is washed out. In such 
investigations, we incorporated the structure of the molecule in the prefactor, and 
kept the same action as in the single-atom case. 

A legitimate question is, however, how the different rescattering scenarios, 
involving one or two centers, contribute to the above-mentioned patterns in NSDI. 
In the following, we will address this issue, incorporating the prefactors derived in |19| 
in the semiclassical action. A closely related procedure has been followed in [71 [9] , 
for high-order harmonic generation. In order to facilitate this assessment and make 
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the two-center interference fringes as clear as possible, we strip the problem to its 
bare bones using a series of assumptions. Firstly, unless otherwise stated, we consider 
very simplified highest occupied molecular orbitals, for which analytic interference 
conditions have been derived [19j. Furthermore, we choose the interaction through 
which the second electron is dislodged as a contact-type interaction at the positions 
of the ions. As it will be discussed subsequently (c.f. Sec. [3]), this type of interaction 
eliminates any bias in the electron momentum distributions, which may be detrimental 
for their interpretatiorlj. Finally, throughout, we will consider the condition for which 
the interference patterns are most pronounced, i.e., parallel alignment. 

This paper is organized as follows. In Sec. [2l we give the expression for the 
transition amplitude related to the process in which the second electron is freed 
by electron-impact ionization, within the SFA. We also briefly recall the two-center 
prefactors employed in [19], and the expression derived in Ref. [19] for the two- 
center NSDI interference conditions (Sec. 12. 2p . Subsequently, in Sec. 12.31 we provide 
the expressions for the modified saddle-point equations, in which different scattering 
processes are taken into account. In Sec. jS] these equations are employed to compute 
differential momentum distributions, which are analyzed in detail. Finally, in Sec. JH 
we summarize the paper and provide its main conclusions. 



2. Transition amplitudes 

2.1. General expressions 

The simplest process responsible for laser-induced nonsequential double ionization 
is electron-impact ionization. Within the Strong-Field Approximation, the 
corresponding transition amplitude is given by 

dt dt' / d3fcy|p,^j}feyfcoexp[i5({P(„)},k,t,t')], (1) 

-OO '^ — oo -^ 
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where the indices in parentheses are related to each electron involved and the term 
{P(n)} in brackets is related to both electron momenta. Eq. ([T]) describes the process 



in which an electron, initially in a bound state 



at a time t' into a Volkov state 



is freed by tunneling ionization 



k(t) ) . Subsequently, this electron propagates in the 
continuum from t' to a later time t. At this time, the electron collides inelastically with 

X Apart from that, in the single-atom case and within the SFA framework, this type of interaction led 
to a better agreement with the existing experiments than more realistic choices, such as a Coulomb- 
type interaction. For a detailed discussion, see Refs. |24II25| . 
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its parent molecule, and a second electron, which is bound in (J)q ), is then released 
through the interaction F12. Thereafter, both electrons are in Volkov states, and their 
final momenta are P(„)(?t. = 1,2). In the above-stated equations, i?o(n)('T' = 1,2) give 
the first and second ionization potentials, and V the atomic potential. The form 
factors ^ and ^ contain all the information about the atomic potential, and the 
interaction by which the second electron is dislodged (see [Ml for details) . One should 
note that the transition amplitude cannot be factorized into one-electron transition 
amplitudes. Physically, this implies a strongly correlated two-electron process. 

Another noteworthy issue is that, within the SFA, T^-o and V/p j^ are gauge 

dependent. In fact, in the length gauge p(„) (t) ~ p(„)+A(r) and k(T) — k-|-A(T)(r = 
t,t'), while in the velocity gauge p(„) (t) — p(„) and k(r) = k. A similar gauge 
dependence has also been reported for high-order harmonic generation [71 [S] and 
above-threshold ionization [18l [11] . For a more general discussion see [22] . 

2.2. Two-center prefactors 

In the specific case of diatomic molecules, we will consider the same simplified model 
as in [in], namely frozen nuclei, the linear combination of atomic orbitals (LCAO) 
approximation, and homonuclear molecules. Under these assumptions, the molecular 
bound-state wave function for each electron reads 



/,(»)^ 



a(")/'». v> /o^ I ^w,(") 



ro'i^(n)) = C^ (/)r(r(„)-R/2) + e0r(r(„)+R/2) , (5) 



where n = 1, 2, e = ±1, and G^ = 1/^2(1 -h eS'(R), with 

S{R) = J [4"^ (r(„) - R/2)] * 0(") (r(„) -I- R/2)dV. (6) 

The positive and negative signs for e correspond to bonding and antibonding orbitals, 
respectively. 

The molecular binding potential, as seen by each electron, is written as 

T/(r(„)) = Fo(r(„) - R/2) + Fo(r(„) + R/2), (7) 

where Vq corresponds to the binding potential of each center in the molecule, which, 
at this stage, are kept general. The above-stated assumptions yield 

V^^ = - (1^ cos[k(t') . R/2]J(k(i')) (8) 

or 



^ko^ = -J^ «in^(*') • R/2]^(k(i')), (9) 

for the bonding and antibonding cases, respectively, with 

J(k(i')) -y'dV(i)exp[zk(t') ■r(i)]yo(r(i))<^^'^(r(i)). (10) 

In the above-stated equations, we have neglected the integrals for which the binding 
potential and the bound-state wave function are localized at different centers in the 
molecule, due to the fact that they are very small for the parameter range of interest. 
If the electron-electron interaction depends only on the distance between the two 
electrons, i.e., V12 — ^(r(i) — '"(2)), the prefactor T^{p, ,}k reads 

C(.,}k = (|^^(p(i) " ^) ^°^[^w • ^m^l^\nt)) (11) 
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or 

<l,}k = (|^^(p(i) - 1^) «"^[^w • ^m^''o\nt)), (12) 

with V{t) — P(i)(i)+p(2)(i) — k(i), for bonding and antibonding orbitals, respectively. 
Thereby, 

^o\'P{t)) = y'dV(2)exp[zP(i) •r(2))]<^^'^(r(2)), (13) 

and 

T/(p(i) - k) = y d^rVir) exp[i(p(i) - k) • r], (14) 

with r = r(i) — r(2). In the velocity and length gauges, the argument in Eqs. (fTTjl . (fT2| 
is given by 'P(t) = p(i) + p(2) - k and 'P(i) = p(i) + p(2) - k + A(i), respectively. 

In terms of the momentum components P(n)||j or p(„)± {n = 1,2), parallel and 
perpendicular to the laser-field polarization, condition ([TT]) or (fT^ may be written as 
cos[(C|| + Ci)-R/2] or sin[(C|| + C±_)R/2], respectively, with 



Cii 

and 



■ 2 



COS 6*, (15) 



C± = P(i)j^ sin6'cosi^ + P(2)± sin0cos(iy9 + a). (16) 

Thereby, the term K{t) is equal to fc ^ A{t) in the length gauge and to k in the velocity 
gauge. Eq. (fT5|) provides well-defined interference fringes, as functions of the parallel 
momenta (p(i)||,P(2)||)- Eq- (HH), on the other hand, has no obvious dependence on 
the alignment angle 6. Indeed, because it depends on the angles (p and a, in the 
momentum plane spanned by the perpendicular momentum components, when one 
integrates over such variables, it main effect is to is blur the interference fringes. In 
this work, we will consider parallel-aligned molecules. This implies that Eq. p^ 
vanishes, and therefore that the interference patterns are sharpest. Explicitly, the 
interference conditions will be given by 

Nn . . 

P{l)\\+P{2)\\^ ^+K{t), (17) 

where N is an integer. For a symmetric combination of atomic orbitals, even or odd N 
gives the interference maxima and the minima, respectively, while in the antisymmetric 
case the situation is reversed. 

2.3. Modified saddle-point equations 

We will now incorporate the structure of the molecule, which is embedded in the 
prefactors ((8|)- p2| . in the semiclassical action. Subsequently, the transition amplitudes 
obtained will be computed employing a uniform saddle-point approximation (c.f. |23j 
for details). For that purpose, we will exponentialize the prefactors Vko and Vjp }k- 
This procedure allows one to single out different rescattering scenarios, involving one 
or two centers. 

This yields the sum 

2 2 

Af = ^^Af,, (18) 
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of the transition amplitudes 

oo t 

Af,„= [ dt [ dt' fd^k^^^\nt))I{k{t'))eMiSj,{{Pin)},Kt,t%(19) 



-OO — oo 



with the modified action 5'jv({p(„)}, k, i, i') defined as 

s.AiPin)}, k, t, t') = 5({p(„)}, k, t, t') + {-iY[r{t)+{-iy+^W)\ • y, (20) 

where the indices j, v relate to the centers in the molecule. In the above-stated 
equation, for processes involving a single center, j + z/ is even, while for scattering 
scenarios in which the first electron leaves at a center Cj and rescatters with a 
center C^,v / j, j + u is odd. The saddle-point equations 9fS'jv({p(„)},k, i,i') = 
5t'5'ji.({p(„)},k,t,t') = and ak5'ji,({p(„)},k,t,i') = are explicitly given by 

^^^^^^'^^' = -i?o(i) + {-ir+^+'dMt') ■ R/2, (21) 

ft 



y^ dr [k + A(r)] + (-l)'^9k { [(-l)-''+''k(t') - k(t)] ■ R/2} = 0, (22) 



and 



n=l 

Eq. (|2ip expresses the conservation of energy at t', with tunneling ionization of the 
first electron. Eq. ([22)) provides the condition for the first electron to return, either to 
the site of its release or to the other ion. Finally, Eq. (P5)) yields the conservation of 
energy at the instant of rescattering. 

One should note that the saddle-point equations (PT|) and (^51) are gauge 
dependent. Specifically, in the length gauge, the tunneling and rescattering conditions 
are given by 

^^^^^^^'^^' = -i?o(i) + (-l)2'^+-'E(i') . R/2 (24) 

and 

respectively. Physically, the additional terms in Eqs. (|24p and (P5|) may be interpreted 
as potential-energy shifts, which sink or increase the ionization potentials. For 
instance, for single-center processes (j = v), such shifts are symmetric, whereas for 
two-center scattering scenarios [j ^ v) they possess the same sign. In the velocity 
gauge, the corresponding equations read 

[k + A(t')]'--2i?o(i), (26) 

and 

^ [P(„)+A(t)p _ [k + A(^)]^ 

n=l 
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respectively. Eq. ([^5)1 and (P^)) are identical to the NSDI saddle-point equations for 
a single atom expressing tunneling ionization at t' and conservation of energy at t. 
Furthermore, in the limit i?o(i) ~^ 0: Eq. P^ may be interpreted as the initial kinetic 
energy of a classical particle, which starts its motion with vanishing drift momentum. 
If written in terms of the electron momentum components parallel and 
perpendicular to the laser-field polarization, the rescattering conditions ((25|) and ([27| 
read 



5 1 5 = 7^ -^^0(2) ■ (.^8 j 



E 

ri=l 

In the length and velocity gauges, i?o(2) = ^o(2) + (^l)''E(t) • R/2 and £'0(2) — £'0(2), 
respectively. In the former case, £'0(2) can be viewed as an effective second ionization 
potential. The above-stated expression gives the equation of a six-dimensional 
hypersphere in momentum space, whose radius corresponds to the momentum region 
for which electron-impact ionization has a classical counterpart. If the kinetic energy 
of the first electron, upon return, is smaller than £9(2) i then the second electron cannot 
be released and this process is classically forbidden. In this case, the corresponding 
transition probability is vanishingly small. 

The return condition, and how it is related to such indices, can be clearly seen in 
Eq. (22). In fact, if j = v, the additional terms vanish, and the condition 

dr [k + A(r)] = 0, (29) 



is obtained. Eq. (P5)) constrains the value of the intermediate momentum k so that 
the electron is leaving and returning to the same center in the molecule. On the other 
hand, if j ^ i^, Eq. ^ reads 



dT[k + A(r)] + (-l)'^+iR = 0. (30) 

If 

The negative and the positive sign corresponds to the situation in which the electron 

starts from the center Ci and rescatters with C2 , or starts from C2 and rescatters with 

Ci, respectively. Both Eq. (P5|) and Eq. (P^ are identical to the return conditions 

obtained from the classical equations of motion of an electron starting its orbit with 

vanishing drift momentum, and propagating under the influence of solely the laser 

field. 

One should note that, within the context of high-order harmonic generation and 
above-threshold ionization, the additional potential-energy shifts present in the length- 
gauge SEA have raised a great deal of controversy. Indeed, there has been considerable 
debate whether such shifts are not an artifact of the strong-field approximation. 
Their existence is mainly related to the fact that the ions in the molecule are at 
a distance ±R/2 from the origin of the coordinate system, which is located at the 
geometric center of the molecule [26l[28]. Furthermore, their presence does not allow 
an immediate connection with the classical equations of motion of an electron in an 
external laser field, as in the velocity-gauge formulation. For the above-stated reasons, 
it was suggested that these additional phase shifts should be removed by dressing the 
electronic bound states [28] . 

In recent HHG computations, however, it has been found that, if the additional 
potential-energy shifts are removed by employing field-dressed states, the two-center 
interference patterns break down [9l [27]. A direct comparison with the numerical 
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solution of the time-dependent Schrodinger equation has shown that this should only 
occur for very large internuclear distances, for which the atoms constituting the 
molecule practically behave as isolated entities [21|- Unfortunately, for NSDI, there 
are neither realistic computations by other means nor enough experiments to settle the 
issue. Therefore, in the following, we will compute electron momentum distributions 
in both length and velocity gauges without favoring any of them. 

3. Electron momentum distributions 

We will approximate the external laser field by a monochromatic wave E(i) — 
So sin ujtex. In this case, the electron momentum distributions, as functions of the 
momentum components (p(i)||,P(2)||) parallel to the laser- field polarization, read 

i^(P(l)|hP(2)||) = yy rf'p(l)±rf'p(2)±|Mfl(p(l),P(2)) + A/L(P(l),P(2))|',(31) 

where the indices R and L denote the right and left peaks, respectively. The amplitude 
Mr is given by Eq. (P), and Ml(p(i),P(2)) = Mi?(-p(i), -p(2)). In Ml (p(i),p(2)), 
the action, the field and prefactors are displaced by half a cycle with respect to Mr, 
i.e., we used the symmetry A(i) = ±A(i ± T/2). 

The above-mentioned symmetry will guarantee that the overall electron 
momentum distributions, obtained either using the prefactor (|5))- p^ or all four 
different transition amplitudes, will be symmetric with regard to the simultaneous 
inversion (p(i)||,P(2)||) — * i~P(i)\\j ~P{2)\\) of the parallel momentum components. 
Specifically for the transition amplitudes M^v, we observe that Mjv(p(i),p(2)) = 
M^j(— P(i), — P(2)). If both start and rescattering times t,t' were real, this would 
lead to symmetric distributions upon (p(i)||,P(2)||) -^ (— P(i)||i ~P(2)||)- Furthermore, 
this would also guarantee that \Mii + Mi2|^ — IM22 + A/2i|^ and |A/i2 + Af22p = 
\Mii + Af2ip. Since, however, the first electron is released through tunneling 
ionization, Im[t'] ^ and this is not necessarily so. This issue is particularly important 
in the length-gauge formulation of the SFA, and will be discussed subsequently. 

In the form factor VkOj we considered Coulomb-type binding potentials Vq and a 
symmetric combination of Is orbitals. We assume that the second electron is dislodged 
by a contact-type interaction T/12 placed at the position of the ions. Explicitly, 

F12 = (5(r(i) - r(2)) [,5(r(2) - R/2) + 5(r(2) + R/2)] , (32) 

and can be viewed as an effective interaction, which in a rough way, accounts for the 
presence of the residual ions. This interaction has been employed in Ref. [19J, within 
a two-center context, and has also been widely used in the single-atom case [Ml I25| . 
In the specific context of this work, it has the advantage of eliminating any additional 
momentum dependence from (p^ {V{t)), whose effect could be hard to disentangle 
from the two-center interference. 

Throughout, we employ a higher driving- field intensity than those reported in 
typical NSDI experiments involving diatomic molecules [IHl IS] ■ This has been done 
with the purpose of extending the region in momentum space for which electron- 
impact ionization exhibits a classical counterpart, i.e., the radius of the hypersphere 
(|28p. In previous work, we have shown that this was necessary in order to make a 
detailed assessment of quantum-interference effects in this context [19] . Furthermore, 
we restrict ourselves to the case of a symmetric combination of Is orbitals. This gives, 
for the specific interaction ((32)) . 

nP(„,}k ~ cos[r{t) ■ R/2]Vioo(0) (33) 
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Figure 1. Contributions from different scattering scenarios to the electron 
momentum distributions as functions of the momentum components (p(i)||,P(2)||) 
parallel to the laser-field polarization. The distributions have been computed in 
the velocity gauge and for a symmetric combination of Is orbitals. The field 
intensity and frequency have been taken as I = 1.5 X 10^*W/cm , and a; = 0.057 
a.u., respectively, and the ionization potentials Eqi = 0.57 a.u. and Eo2 = 0.98 
a.u. correspond to N2 at the equilibrium internuclear distance _R = 2.068 a.u. 
The upper panels display the contributions from topologically similar scattering 
scenarios, involving only one or two centers, i.e., the transition probabilities 
|Mli -I-M22P and |Afi2 + A/21 P [panels (a) and (b), respectively]. In panels (c) 
and (d), we display the contributions from the processes starting and ending at 
center Ci, respectively (transition probabilities |A/ii + A/i2|^ and |A'/ii + M2i|^, 
respectively). Panels (e) and (f) depict the contributions |Af2i + Af22p ^md 
lAfi2 + Af22p from those starting and ending at C2, respectively. The sum 
Afii + M12 + Af2i + Af22p of all contributions are provided in panel (g). For 
comparison, panel (h) has been computed using the symmetric prefactors ^ and 
mil l and single-atom saddle-point equations. 



Laser-induced nonsequential double ionization in diatomic molecules 10 

Without loss of generality, however, our studies, as well as the conclusions of this 
work, may be extended to the antisymmetric case, or to more complex orbitals. In 
particular, if the interaction ([5^ is taken, only the contributions from s states would 
be non- vanishing in the form factor V{p, ,}k- 

In Fig. [U we display the electron momentum distributions in the velocity gauge 
and for symmetric orbitals, considering specific scattering processes. For comparison, 
in Fig. [11(h), we are presenting the distributions obtained employing the prefactors 
([8]) and (fTTIl and solving single-atom saddle point equations. This approach has been 
considered in |T^, and incorporates all the structure of the molecule in the prefactors. 
This distribution exhibits fringes in agreement with the interference condition (|17p . 

In Fig. [TJ(a), only the contributions from the processes in which the first electron 
leaves and returns to the same center have been taken into account. These processes 
are related to the transition amplitudes Mjj{j = 1, 2). In this case, fringes parallel to 
the anti-diagonal p(i) 1 1 = —p{2)\\ are also present. Their position, however, disagrees 
with condition (fT7)) . A closer inspection, however, shows that the terms Mn and M22 
can be grouped as cos[(p(i) -I- P(2)) ■ R-/2]. This gives another interference condition in 
terms of the parallel momentum components, namely P(i)|| + P{2)\\ — N-k / R. Maxima 
and minima are present for even and odd N, respectively. We have verified that 
all maxima and minima in the figure are approximately given by the above-stated 
expression, and correspond to the integers < TV < 4. 

If only the two-center processes are taken [Fig. [U(b)], the transition amplitudes 
M21 and M12 can be grouped as cos[(p(i)-|-p(2)— 2k)-R/2]. This gives fringes following 
P(i)\\ +P(2)|| — '2,k = Ntt/R, which, once more, differ from the global interference 
condition (fT7|) . A rough analytic estimate of the position of the fringes can be made, 
by considering that the first electron leaves at peak field and returns at a field crossing. 
In this estimate, instead of using the modified return conditions ([30)) . we considered 
the saddle-point equation (|29p. which states that the electron returns to the site of its 
release. In the present context, this expression for k may also be viewed as an average 
value between those given by the two-center return conditions given in ([30]) . For the 
specific parameters in this work, this yields P(i)\\ +P(2)|| — (1-45A'' — 0.849)^yUp. This 
estimate is in good agreement with the maxima and the minima displayed in Fig.[T](b). 

Interestingly, the distributions obtained from the transition probability \Mjj + 
Mjvp, with j = 1,2 and j ^ z/, agree with the overall interference condition pT]) . 
This is shown in Figs. [T](c) and[T](e), and is due to the fact that one may rewrite the 
sum of such terms as exp[±ik • R/2] cos[(p(i) -I- p(2) — k) • R/2] . Apart from an overall 
phase factor, this is the same interference condition as if all transition amplitudes are 
taken. This latter case is shown in Fig. [TJ(g), for comparison. Physically, the term 
Mjj corresponds to the orbits in which the first electron is freed at and subsequently 
rescatters off a specific center Cj, while Mj^, gives the process in which it reaches 
the continuum at Cj and returns to the other center (7,^. This suggests that these 
processes are the most relevant in determining the interference patterns from NSDI in 
a molecule. 

On the other hand, if we consider scattering scenarios ending at the same center, 
the electron-momentum distributions resemble very much those obtained for the 
single-atom case, namely isotropic distributions centered nearp(i)|| = P(2)\\ = ±2-^/^/^^. 
This holds regardless of whether one or two centers are involved, and can be seen 
in Figs. [11(d) and [T](f). In each panel, both one and two center processes are 
considered. However, interference fringes are absent. This is consistent with the 
double-slit physical picture, which relates the existence of the interference maxima 
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Figure 2. Contributions from different scattering scenarios to the electron 
momentum distributions as functions of the momentum components (P(i)|| ,P(2)||) 
parallel to the laser-field polarization. The distributions have been computed in 
the length gauge and for a symmetric combination of Is orbitals. The remaining 
parameters are the same as in the previous figure. The upper panels display the 
contributions from topologically similar scattering scenarios, involving only one or 
two centers, i.e., the transition probabilities IM11-I-M22P and \Mi2+M2i\^ [panels 
(a) and (b), respectively]. In panels (c) and (d), we display the contributions 
from the processes starting and ending at center Ci, respectively (transition 
probabilities |Mii -I-M12P and |Mii -I-M21P, respectively). Panels (e) and (f) 
depict the contributions | Af2i + M22 P and |Mi2 + M22 P from those starting and 
ending at C2, respectively. The sum \M11-\-M12 -\-M21 -I-M22 P of all contributions 
are provided in panel (g). For comparison, panel (h) has been computed using 
the prefactors l[8) and JTTJ and the single-atom saddle-point equations. 
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and minima to photoelectron emission in spatially separated centers. This can also be 
seen by rewriting these contributions as exp[±z(p(i) + p(2) — k ■ R/2)]Vko- The form 
factor V]t:o does not lead to well-defined interference patterns, so that the distributions 
are very similar to those obtained in the single-atom case (for discussions see [TOll^ V 
Finally, we observe that Figs. [T](d) and[TJ(f) are practically identical. This means 
that the condition |Mii -t- M21P = |Af22 + M12P holdfU. Intuitively, this is expected, 
as the transition amplitudes Mji{j — 1,2) are the complex conjugates of Mj2- The 
same holds for Figs. [11(c) and[TJ(e), which expresses the fact that \Mii + Mi2\^ — 
IM22+M21I2. 

The length-gauge counterparts of the distributions discussed above are presented 
in Fig. [2l In this case, we expect a different interference condition, according to 
Eq. (fTT]) . Furthermore, an inspection of the saddle-point equations shows that there 
exist additional potential-energy terms which lower or increase the barrier through 
which the first electron tunnels out. These energy shifts will influence the electron- 
momentum distributions, as we will discuss subsequently. 

In the upper panels, we consider topologically similar scattering scenarios, 
involving only one or two centers [Figs. [2] (a) and[2](b), respectively]. In contrast 
to the velocity-gauge situation, these distributions exhibit at most a slight distortion, 
as compared to the single-atom case. Sharp interference fringes, however, are absent. 
This is a consequence of the above-mentioned potential-energy shifts. Such shifts alter 
the potential-energy barriers in such a way that is much more probable for the first 
electron to tunnel out from a specific center of the molecule. Hence, the processes 
starting from a center Cj in the molecule will be far more prominent than those 
starting from the other center Cu{v ^ j) and there will be no sharp fringes. For 
the parameters employed in this figure, we estimate a difference of roughly one order 
of magnitude between each transition amplitude. A similar effect was present for 
high-order harmonic generation and has been discussed in detail in Ref. [9] . 

Also in the length gauge, the transition amplitudes Mjj + Mj^, related the 
processes starting at the same center and ending at different centers lead to the same 
interference condition as the overall prefactors ^ and HI]), apart from a phase factor 
exp[±j^], with ^ = [k -I- A{t')] ■ R/2. There are, however, some major quantitative 
differences, due to the fact that Im[i'] ^ 0. This leads to a non-vanishing imaginary 
part in ^, which, depending on the field, will enhance the yield for a center Cj and 
suppress it for the other center Ci, in the molecule. As the field reverses its sign, the 
contributions from Cj will be suppressed and those from the other center enhanced. 
This leads to strongly asymmetric electron momentum distributions, and may be 
interpreted as another consequence of the additional potential-energy shifts, which 
are present only in the length-gauge formulation of the SFA. 

The above-stated features are explicitly shown in Figs. [21(c) and[2](e), starting 
from Ci and C2, respectively. In the former case [Figs. [2] (c)], the yield is strongly 
suppressed in the region of positive parallel momentum. This is due to the fact that 
the additional potential energy shifts suppresses tunneling ionization for all the orbits 
starting from Ci. As the field reverses its sign, i.e., when t' ^ t' + r/2, the barrier 
will be sunk by these shifts. This explains the much larger yield in the negative 
momentum region. Fig. [21(e), which is the mirror image of Figs. [2] (c) with respect 

§ We have verified, by a series of numerical computations, that this condition only holds in practice. 
This is a consequence of the fact that the start time t' exhibits a non- vanishing imaginary time, and 
that, according to the saddle-point Eqs. (29), (30), the intermediate momentum depends on this 
variable. The relative differences in the yields, however, is less than 0.2%. 
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to (pi||,P2||) ^ {^Pi\\j^P2\\)^ exhibits the consequences of the opposite effect: for 
the time t', the barrier is sunk for C2. This leads to enhanced distributions in the 
positive momentum regions. After half a cycle, as the field changes its signs, the 
barrier will be lifted and and the contributions in the negative momentum region 
suppressed. A similar effect has been observed for high-order harmonic generation 
[SI Hn], In particular, in |10j . apart from modified saddle-point equations, we have 
employed effective prefactors and single-atom saddle-point equations in order to mimic 
the interference between different processes. Also in the latter case, we observed 
different orders of magnitude for processes starting at different centers. 

Therefore, in contrast to the velocity- gauge situation, the processes starting at 
a particular center of the molecule do not lead to the same distributions as if all 
four terms in Eq. (|20|) are taken. These latter results are displayed in Fig. [21(g). 
For comparison, in Fig. [21(h), we also provide the outcome of the computations 
using single-atom saddle-point equations and the prefactor pTjl and ([5]). Finally, 
if the processes starting at different centers, but ending at the same center lead to 
distributions qualitatively similar than those observed in the single-atom case [see 
Fig. [21(d) and (f)]. Similarly to the velocity-gauge situation, the terms can be 
combined in exp[±i(p(i) -I- p(2) — (k — A{t)) ■ R/2)] cos[k -I- A(t')]. Even if the extra 
dependence on the vector potential influences the distributions, qualitatively, similar 
conclusions hold. However, there is a slight asymmetry with respect to the reflection 
iP(i)\\jP2\\) ~^ (~P(i)||i ~P(2)||) due to the fact that the times i,i' are complex. Once 
more, the contributions of the processes ending at a particular center Cj are the mirror 
image those ending at the other center Ci,. 

In the following, we explicitly show that the above-mentioned asymmetries are 
caused by Im[i], and Im[i']. For that purpose, in Fig. 3, we plot the counterparts 
of Fig. 2(c)-2(f) employing only the real parts of such times. In this case, the 
contributions of the orbits starting at the same center are symmetric, i.e., the condition 
|Afii + M12P = lA/22 + M21P holds [see Figs. 3.(a) and 3.(c)]. The contributions 
ending at the same center, depicted in Figs. 3.(b) and 3.(d), also exhibit the symmetry 
|Mii -|- M21P = IM22 + Mi2p. In these computations, instead of employing modified 
saddle-point equations, we grouped the separate transition amplitudes in the effective 
prefactors 

VSf'^\{p^,,j},k) ~ exp[z(-l)^"[k + A(<')] • R/2] (34) 

X cos[(p(i) + p(2) - k + A(<)) • R/2] 

and 

^cf ({P(n)},k) ~ exp[*(-l)^"(p(i) +p(2) -k-fA(t) •R/2)] (35) 

X cos[k + A(t')] 

mimicking the processes starting and ending at a center Cj , respectively, and employed 
single-atom saddle-point equations (see Ref. [TD] for details on this procedure). The 
imaginary parts of i, i' were removed in the prefactors only. 

4. Conclusions 

We have made a detailed analysis of the different scattering scenarios in laser-induced 
nonsequential double ionization of diatomic molecules, within the framework of the 
strong-field approximation. We assumed that the second electron is dislodged by 
electron-impact ionization and employed saddle-point methods. The semiclassical 
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Figure 3. Transition probabilities starting and ending at the same centers, 
considering only the real parts of t, t', as functions of the momentum components 
(P(i)||iP(2)||) parallel to the laser-field polarization, for the same parameters as in 
Fig. 2. Panels (a), and (b) exhibit the contributions from the processes starting 
and ending at center Ci (transition probabilities |Mii -I-M12P and |Mii -I-M21P, 
respectively), while panels (c) and (d) depict the contributions IM21 -I-M22P and 
\Mi2 + Af22p from those starting and ending at C'2, respectively. 



action has been modified in order to account for four different physical processes. 
In two of them, the first electron is released and rescatters off the same center 
Cj{j = 1,2). In the remaining processes, it leaves from a center Cj and rescatters 
with a different center Cy^ v ^ j. 

We placed particular emphasis on the quantum interference between such 
processes, and on how it affects the differential electron momentum distributions, 
as functions of the momentum components p(„)||(n = 1,2) parallel to the laser- 
field polarization. For that purpose, we considered parallel-aligned molecules, for 
which the interference patterns are sharpest [TS]. We found that the contributions of 
the processes in which the first electron starts at the same center Cj, regardless of 
whether it returns to the site of its release or to the other center C„ , yield interference 
patterns which fulfill the overall interference condition (flT)) . There exist, however, 
differences for both velocity and length-gauge formulations of the SFA. Indeed, while 
in the velocity gauge the contributions of the above-stated orbits leads to the same 
distributions as in the overall case, this affirmative cannot be made in the length-gauge 
situation. In this latter case, there are additional potential energy shifts, which sink 
or increase the potential barrier through which the first electron tunnels. This favors 
the contributions of a particular center in the molecule, and suppresses those from 
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the other center. As the electric field changes sign, so do the shifts. This leads to 
asymmetries in the electron-momentum distributions. 

The quantum interference between other types of processes leads either to quite 
different interference patterns, or to momentum distributions that strongly resemble 
those obtained for a single atom. In particular, a resemblance to single-atom 
distributions is observed if only the processes for which the first electron rescatters 
at the same center Cj are considered. The fact that interference fringes are absent is 
in agreement with the double-slit physical picture, as there is no electron emission at 
spatially separated centers for such processes. 

Sharp interference fringes, on the other hand, are observed if we consider 
topologically similar scattering scenarios, involving either one or two centers, and 
the velocity gauge. This is expected since, in this case, the first electron rescatters 
at spatially separated centers. It is worth noticing, however, that the interference 
patterns, and thus the electron-momentum distributions, look quite different from 
those obtained if all contributions are taken. Physically, this indicates that we can 
not single out the interference between such processes as being the most relevant. 
In the length gauge, in principle, there are also interference patterns. In practice, 
however, we only see slight distortions, as compared to the single-atom case. This 
is due to the additional potential-energy shifts, which sink or increase the potential 
barrier through which the first electron must tunnel in order to reach the continuum. 
Hence, the contributions starting at a center Cj in which the barrier has been sunk 
will be far more prominent than those from a center C^, in which it has been raised. 

The above-stated issues contribute to the present controversy involving the above- 
mentioned potential-energy shifts [H El [Til UHl EH EH EH]- Indeed, it is not clear 
whether they are only an artifact of the strong-field approximation, or whether they 
possess some physical meaning. For instance, in jHEZ] it has been shown that, if such 
shifts are not present, there exists a breakdown in the double-slit interference patterns 
for high-order harmonic generation. Specifically in [27], it was argued that, for small 
internuclear distances, there is no physical justification for dressing the electronic 
bound states with a linear phase factor. In fact, according to [57j a quadratic phase 
factor would be the correct form of dressing, and would physically correspond to the 
quadratic Stark effect. On the other hand, in [7], it was shown that, for internuclear 
distances of the order of the electron excursion amplitude, the additional potential- 
energy shifts led to an unphysical increase in the high-harmonic intensities and cutoff 
energies. Furthermore, in [28j it was demonstrated that a linear dressing in the 
electronic bound states is necessary in order cancel out the extra potential-energy 
shifts and to restore, among others, translation invariance. In this specific work, one 
argument against the existence of the potential-energy shifts would be that, for a 
homonuclear molecule, one would intuitively expect the symmetries jMn -I- M21P = 
\M22 + M12P and jMn + M12P = IM22 + M21P to hold. Due to the fact that, in the 
length gauge, the imaginary part of t' is shifted in different ways for different centers 
in molecule, we have seen that, instead, the distributions starting at a specific center 
Cj are a reflection upon (pi||,P2||) ^ i^Pi\\i ^P2\\) of those starting at the other center 

Finally, we would like to comment on the similarities and differences between 
the present results and those reported in our previous work [HI [IH] , on the double-slit 
quantum interference for HHG in diatomic molecules. Therein, even though there 
was strong evidence that the interference maxima and minima in the spectra were 
due to the transition probability \Mjj + Mjv\^, related to the interference of the 
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orbits starting at the same center and finishing at different centers, there were several 
ambiguities. Firstly, this information could only be extracted in the length gauge, 
due to limitations of the SFA. Within this approach, there is a breakdown in the 
interference patterns if the velocity gauge is employed to compute the high-harmonic 
yield. Furthermore, similarly to the present case, the length-gauge formulation of the 
SFA exhibits additional potential-energy shifts which make the ionization probability 
far more probable from one center than from the other. This holds even for very small 
internuclear distances [lOj . and could be the reason why the interference patterns 
would not be present also in the transition probabilities iMjj+AI^jyl'^, or |Mjv + Af^jp, 
involving one or two-center orbits, respectively. By providing an additional pathway 
for the electron to reach the continuum, using attosecond pulses, we could verify that 
these processes led to maxima and minima if the potential-energy shifts could be 
overcome. Since, however, we were modifying the physics of the problem, we could 
not reach a definitive conclusion. 

Fortunately, many of these ambiguities are not present for laser-induced 
nonsequential double ionization. This is possibly due to the fact that this 
phenomenon intrinsically involves electron-electron correlation, whereas high-order 
harmonic generation may be modeled within the single-active electron approximation. 
Indeed, although the distributions computed in this work employing the topologically 
similar scattering scenarios may, in some cases, exhibit interference patterns, these 
do not agree with the overall interference conditions. Furthermore, in contrast to the 
high-order harmonic case, there is no breakdown of the interference patterns in the 
velocity-gauge situation. This allows one to assess the influence of the topologically 
similar scenarios on the NSDI electron momentum distributions in the absence of the 
above-mentioned potential energy shifts. Finally, it is rather interesting that, in the 
length gauge, the influence of the potential energy shifts can be directly mapped into 
the electron momentum distributions for the processes starting at the same center. 
For instance, if the electron started at a center Ci, for the parameters used in this 
work, there is a suppression of the yield in the positive momentum region and an 
enhancement for positive parallel momentum, whereas if it has started at C2 the 
situation will be reversed, i.e., there will be an enhancement in the region of positive 
parallel momenta and a suppression forpj^jn < {n = 1,2). This means that different 
potential energy shifts lead to an enhancement or a suppression of the distributions 
in distinct momentum regions. 

Hence, if these processes could be isolated in a realistic scenario, one could 
in principle determine from which side of the molecule the first electron reached 
the continuum. Furthermore, if no asymmetry is observed, the existence of the 
potential-energy shifts and the validity of the length-gauge SFA electron momentum 
distributions could be ruled out in the context of diatomic molecules. For that purpose, 
it would be necessary to suppress all the the orbits starting from one center of the 
molecule in the overall distributions . This could be achieved if, for instance, the 
system is prepared in such a way that the initial state of the first electron is localized 
in one center of the molecule, or if the external laser field is taken as a few-cycle pulse 
with an adequate choice of carrier-envelope phase. 
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